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We report on the experimental observation of a novel interorbital Feshbach resonance in ultracold 
^^^Yb atoms. This opens up the possibility of tuning the interactions between the ^So and ®Po 
metastable state, both possessing zero total electronic angular momentum. The resonance is ob¬ 
served at experimentally accessible magnetic field strengths and occurs universally for all hyperfine 
state combinations. We characterize the resonance in the bulk via interorbital cross thermalization 
as well as in a three-dimensional lattice using high-resolution clock-line spectroscopy. Our measure¬ 
ments are well described by a generalized two-channel model of the orbital-exchange interactions. 

PACS numbers: 34.50.Cx, 75.10.Dg, 67.85.Lm 


Magnetic Feshbach resonances have become an indis¬ 
pensable tool in the study of ultracold quantum gases, 
enabling the tuning of interaction strengths over a wide 
parameter range [1]. This tunability has given rise to an 
impressive set of experimental results, including the re¬ 
alization of the BEC-BCS crossover in degenerate Fermi 
gases [2-5] as well as the discovery of novel few-body 
phenomena such as Efimov trimers [6]. 

Whereas the majority of alkali atomic species feature 
Feshbach resonances, this is not the case for alkaline- 
earth-like atom states without electronic angular momen¬ 
tum. In this Letter, we report on the observation of a 
recently predicted novel type of Feshbach resonance [7], 
enabling the tuning of interorbital interactions based on 
the Zeeman shift of different nuclear spin states of the 
atoms. 

In a Feshbach resonance, a bound state of the ener¬ 
getically inaccessible molecular potential (closed chan¬ 
nel) couples to the scattering continuum of the energet¬ 
ically allowed entrance channel. This coupling dramat¬ 
ically affects the atomic scattering properties whenever 
the bound state energy is close to the open channel en¬ 
ergy. In magnetic Feshbach resonances, an external mag¬ 
netic field can tune these states into and out of reso¬ 
nance, as they have different magnetic moments such as 
in the prototypical example of singlet and triplet scatter¬ 
ing channels of two alkali atoms [8-10]. 

In contrast, alkaline-earth-type atoms with two valence 
electrons such as ytterbium and strontium possess no 
electronic angular momentum in the atomic ground state, 
either for the electronic spin singlet ^Sq (denoted \g)) or 
for the triplet ^Pq (|e)). This and the associated suppres¬ 
sion of hyperfine coupling make these atoms superb can¬ 
didates for atomic clocks [11] and for the investigation of 
new exotic many-body states with extended SU(A'^) sym¬ 
metry [12]. However, this also implies that no magnetic 
Feshbach resonances are expected within the two orbital 
states. 



FIG. 1. Illustration of the magnetic field dependence of the 
scattering potentials in the two-channel orbital Feshbach res¬ 
onance. Blue solid lines represent the scattering potentials 
at zero magnetic field. The electronic triplet (nuclear spin 
singlet) supports a bound state with energy es. A finite mag¬ 
netic field induces a deformation of the scattering potential 
for large atomic separation (yellow solid lines), proportional 
to the differential magnetic moment. The inset illustrates the 
origin of the differential magnetic moment due to different 
Lande g-factors of the ^So (bright symbols) and ^Po (dark 
symbols) states. 

Instead, optical Feshbach resonances [13], based on the 
coupling to a bound molecular state via a one- or two- 
photon process have been investigated [14, 15]. Because 
of their limited tunability and lifetime, such optical Fes¬ 
hbach resonances have been difficult to exploit in exper¬ 
iments. Until now, magnetic Feshbach resonances in yt¬ 
terbium have only been observed when atoms are specifi¬ 
cally prepared in the ^P 2 state possessing electronic mag¬ 
netic moment and hyperfine coupling [16]. However, re¬ 
cent measurements [17, 18] of the scattering properties 
of a ^So - ^Pq atom pair in two different nuclear spin 
states suggest the existence of a shallow molecular bound 
state, leading to the prediction of a magnetically acces¬ 
sible scattering resonance [7]. 

The scattering resonance is described by interaction 
channels possessing both orbital and nuclear degrees 
of freedom in ^^^Yb (/ = 5/2). As in Scazza et al. 
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[17], the interaction can be expressed in the basis con¬ 
sisting of symmetric and antisymmetric superpositions 
|i) = 5 (l^e) i \^9)) (Iti) T lit)) of the electronic or¬ 
bitals (|e),|( 7 )) and the two nuclear spin states (|t))lt)) 
with m\;,,rn'p G —Associated to |±) are the 
orbital singlet scattering length a~g = 219.5(29) oq [17] 
and the very large triplet scattering length a+ > 2000 oq 
[17, 18], with oq denoting the Bohr radius. The two or¬ 
bitals exhibit different Lande g-factors, due to a weak 
hyperfine coupling of the ^Pq with the ^Pi state, giving 
rise to a differential magnetic moment Apt — (ggTnjp — 
ggm^F)fiB, where gs is the Bohr magneton. The pres¬ 
ence of a magnetic field will therefore mix the singlet 
and triplet states, introducing a coupling between the 
orbital and spin degrees of freedom. In this case, the non¬ 
interacting system has the eigenbasis |o) = ([gtiei) ~ 
|e i; g t))/V2 and |c) = (|e t; 5 i) - Is i; e t))/v^- 

Two atoms entering the potential in the open channel 
eigenstate |o) couple to |c) through the orbital-mixing 
interaction term. At short distances the molecular inter¬ 
action potentials dominate and are independent of the 
magnetic field due to the symmetric nature of the |±) 
states. Nevertheless, the entrance channel |o) can be 
brought into resonance with bound states of the interor¬ 
bital molecular potentials by shifting the entrance energy 
using the differential Zeeman shift AfiB as illustrated in 
Fig. 1. The resonance occurs when AfiB comes close 
to the binding energy of the bound state in the closed 
channel Cb ~ —h^Kma^), with Uc = {a^g -\- a~g)l2. 

To determine the magnetic field dependence of the 
scattering length, we perform cross-dimensional thermal- 
ization measurements in a 3D harmonic trap [19]. A mix¬ 
ture of |e ),) and [g f) is prepared in an out-of-equilibrium 
state and rethermalizes with the mean rate of elastic colli¬ 
sions, which can be derived from Enskog’s equation [20]. 
In the limit of low-energy scattering, the collision rate 
(and therefore the thermalization rate) is proportional 
to the elastic scattering cross section a. 

Our experiments begin by creating a degenerate Fermi 
gas via evaporative cooling in a crossed dipole trap [17]. 
The desired spin mixture is prepared by optical pumping, 
resulting in a two-component \g 1) \g f) gas with typical 
temperatures T ~ 0.2 Tp and a total atom number of 
= 3 X 10"^ per spin state, where Tp is the Fermi tem¬ 
perature. 

In order to populate the ^Pq state, the atoms are 
first loaded into a one-dimensional, state-independent, 
i.e., magic-wavelength Am = 759.3 nm lattice [21] in the 
Lamb-Dicke regime. This lattice is coaligned with the 
579 nm excitation beam used to apply a 7r-pulse on the 
\g ),) — > |e ),) optical clock transition at an applied mag¬ 
netic field of i3exc = 1200 G. Next, the atoms are released 
adiabatically into a magic-wavelength dipole trap (see 
supplementary information [22]). The sample is then 
heated along the z-direction by repeated Bragg pulses 



FIG. 2. Cross-thermalization rate as a function of the mag¬ 
netic field B for one- and two-orbital mixtures. Circles and 
diamonds mark the thermalization rate of the je J,) \g t) and 
the \g 4.) \g t) mixture, respectively with Amp = 5 (mp = 
— 5/2, -1-5/2). Square symbols show the values for a je 4,) [^ t) 
mixture with Amp = 1 {mp — —5/2,—3/2) rescaled to 
B —>■ ^AmpB. Error bars denote the la uncertainty of the 
fit to the cloud aspect ratio. The inset shows a conversion of 
the thermalization rate to the scattering length Uerstj based 
on the reference provided by the ground state Ogg = 199.4 ao 
[25] (dashed line). In the inset, the offset due to the residual 
ground state thermalization rate rgg,res = 2.17s“^ has been 
subtracted (see main text). 


using a standing wave along this axis, and the subse¬ 
quent thermalization of excitation into the orthogonal 
directions is observed by measuring the cloud aspect ra¬ 
tio 0(t) after 13 ms of time of flight. We determine the 
thermalization rate T^g by a single exponential fit to the 
time dependence 9{t) = 1-1- aexp(—Fegt) [24]. 

We observe a change in the thermalization rate over 2 
orders of magnitude as a function of the magnetic field 
as shown in Fig. 2. In contrast, for a spin mixture in the 
ground-state orbital |g), we find the thermalization rate 
Tgg to be independent of the magnetic field. The result¬ 
ing B-field dependence of Feg exhibits the characteris¬ 
tic shape of a Feshbach resonance with a peak position 
Bq = 55(8) G, and a zero crossing at Ba = 417(7) G, 
both determined by quadratic fits within ±15 G (±40 G) 
regions around the resonance (zero-crossing) position, re¬ 
spectively. The excitation process has an efficiency of 
approximately 90 % and the thermalization rate in the 
zero-crossing regime of Uefgt therefore bounded from 
below by the thermalization of the \g f) atoms with the 
residual \g ),) atoms. 

A unique property of this new type of Feshbach res- 
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FIG. 3. Atom loss in a |e f,) |g t) mixture held at a magnetic 
field B near the orbital Feshbach resonance, (a) Number of 
atoms in the state \g f) with mp = +5/2 (circles) and the 
residual \g /) fraction with mp = —5/2 (diamonds). All data 
points represent averages of at least eight individual measure¬ 
ments and points less than 2.5 G apart are binned to reduce 
visual clutter. Error bars indicating the standard error of the 
mean are smaller than the marker size, (b) Relative popu¬ 
lation of |3t)i i-e-i number of atoms in Ig/) normalized to 
Ngi^{B) = Ng — Ngi{B), whoro Ng is the ground-state atom 
number without losses, averaged for fields B > 120 G, and 
Ngi{B) is the residual atom number in \g /). The inset shows 
the time trace of the total atom number decay close to the 
Feshbach resonance. 


onance is that it involves only two channels, which ad¬ 
ditionally are fully SU(A^) symmetric in the absence of 
a magnetic field, making the coupling universal with re¬ 
spect to the choice of mp states. We demonstrate this 
by repeating the measurement with another spin mix¬ 
ture with a different value of the differential magnetic 
moment A/i. When A/i is taken into account by rescal¬ 
ing the magnetic field axis correspondingly, the data sets 
collapse onto a single curve without further adjustments 
(see Fig. 2), demonstrating the universal behavior with 
respect to different mp-staie combinations. 

In order to characterize the loss channel of the reso¬ 
nance, we use the same preparation protocol as for the 
cross thermalization, omitting the heating procedure. Af¬ 
ter a fixed hold time of 150 ms, the remaining ground- 


state atoms are imaged spin selectively using an optical 
Stern-Gerlach technique [26]. As shown in Fig. 3a, when 
reducing the holding field we observe a magnetic field- 
dependent loss of the \g f) population starting close to 
the position of the resonance. However, the minimum of 
the loss feature does not occur until the field approaches 
B = 0. Below 50 G, we observe a significant repopulation 
of the \g I) state, since orbital exchange becomes energet¬ 
ically favorable at low fields [17, 18]. Accounting for this 
exchange, the loss feature still remains strongly shifted 
towards lower magnetic fields compared to the resonance 
position in the elastic channel (Fig. 3b). This is in fact 
similar to observations in other mixtures of fermions [27- 
30], but the shift is possibly enhanced due to the intrin¬ 
sically large size of the shallow bound state. 

At the resonance position, we extract a lifetime of 
Ti/e = 386(9) ms at a temperature T/Tp ~ 0.3 and an 
initial peak density of uq ~ 5 x 10^^ atoms/cm^. As¬ 
suming three-body decay for the je4,)jgt) mixture, we 
obtain a loss rate coefficient of = 7.5 x 10“^^cm®/s. 
The observed loss is comparable to previously reported 
values in alkali Fermi gases, e.g., mixtures [29]. This 
realization of a stable strongly interacting Fermi gas ap¬ 
pears favorable for exploring novel superfluidity phenom¬ 
ena with this Feshbach resonance in Yb [7]. 

In a third experiment, we directly probe the two- 
particle interaction on individual sites of a deep, 3D 
magic-wavelength lattice, as described in [17]. This of¬ 
fers the advantage of determining the interaction shift of 
the atom pair with high resolution using clock-line spec¬ 
troscopy. Furthermore, in such a setting the two-particle 
problem can be separated from possible many-body ef¬ 
fects and allows for precise determination of the interac¬ 
tion energies. The energies of the two-particle states that 
can be excited on the open channel transition at varying 
magnetic fields are shown in Fig. 4. The energy is given 
relative to that of two non-interacting spatially separated 
atoms, Co- The data points correspond to resonance po¬ 
sitions obtained in an mp = ±5/2 spin mixture for an 
isotropic lattice depth of V = 29.8i?r, where is the 
recoil energy of the lattice light. 

The lowest observed branch (circles) has negative in¬ 
teraction energy for B < 400 G. For low fields, this state 
corresponds to the molecular branch of j±). Extrapolat¬ 
ing the measurements to i? = 0, we find an on-site bound 
state energy of E^/h = 32(2) kHz. 

The next higher-lying energy branch was used in [17] to 
determine the value of a~g = 219.5(29) oq. We verify with 
spin-sensitive measurements that this state evolves from 
the antisymmetric superposition state j—) to the (band- 
excited) open-channel state jo) with increasing magnetic 
fields. Using spin-exchange oscillations at low fields, an 
indirect measurement of a/g = 3300(300) oq was also per¬ 
formed [18]. 

To model the on-site le)-]^) interaction, we general¬ 
ize the two-channel description of Ref. [7] to the case 
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FIG. 4. Magnetic field dependence of the \e)-\g) atom pair in¬ 
teraction energy on a lattice site. The values are given relative 
to the energy of two non interacting atoms. The assignment of 
the spectroscopy resonances (points) is based on the observed 
transition strengths, which decrease for higher lying energy 
states. Solid lines are solutions of the two-particle problem. 
Light blue bands indicate the range of variations of the theo¬ 
retical model spanned by varying atg and a^g by 10 %. (Inset) 
Schematic diagram of the excitation geometry. 

of isotropic confinement and energy-dependent collisions. 
In the absence of a magnetic field, the interactions 
in the singlet and triplet channels decouple. Conse¬ 
quently, the energy shift E is determined by the condition 
■v/ZT+^cot [5fg{E + eo)] -f n(i?) = 0, which depends on 
the singlet and triplet scattering phase shifts, and 
the renormalised pair propagator in the lattice site, n(Zl) 
[31, 32]. However, in general we need to take into account 
the coupling introduced by the magnetic field, and this 
converts the problem into a matrix equation. Thus, the 
condition for the interaction energy becomes 

det[T-i(Z;-A/rH/2-beo) + n(Z;)] = 0, (1) 

with n(i?) = n(Zl) jo) (o| -I- Y\{E — AfiB) |c) (c| and 
t-HE) = ^/E{cotd-iE)\-){-\+S+iE)\+){+\). For 
details of the derivation, see [22]. Note that our model 
goes beyond the usual solution of the two-body problem 
in a harmonic trap [31] and accounts for how both open 
and closed channels are modified under confinement. We 
apply a low-energy expansion of the phase shifts up to 
and including effective range. The magnetic field depen¬ 
dence of all measured interaction shifts is very well repro¬ 
duced by our model (solid lines in Fig. 4). A best fit of the 
model with a+g as a free parameter yields a+g = 1878 ao. 
The corresponding effective range = 216 ao is calcu¬ 
lated analytically [33, 34], assuming a long-range van der 
Waals-potential with Cq = 2561 a.u. [35]. The fit uncer¬ 
tainty on a)!g is 37 Og, but we expect that the uncertainty 
from model simplifications such as the lattice and scat¬ 
tering phase expansions are comparable or larger than 


this. To indicate the variability, we plot a range corre¬ 
sponding to ±10% variation of both scattering lengths as 
shaded areas in Fig. 4. The value for a~g = 219.7± 2.2 oq 
obtained applying the model to the B — 0 data from [17] 
is consistent with the one reported there, and has been 
used in this work with effective range = 126 ap. 

The spectroscopy results also enable us to derive a 
Feshbach resonance position Bq = G and the zero¬ 
crossing Ba = 327g2 G for the bulk at T = 290 nK. 
Both values are in good agreement with the results ob¬ 
tained from the independent cross-thermalization mea¬ 
surements. 

In conclusion, we have observed a new type of scatter¬ 
ing resonance between different atomic orbitals of ^^^Yb 
arising from strong interorbital spin-exchange interac¬ 
tions. Because of the SU(A)-symmetric nature of the ex¬ 
change interaction [17], the resonance occurs universally 
for any spin combination. We have precisely character¬ 
ized this Feshbach resonance using a new model of the 
on-site atom pair interaction energy shift that incorpo¬ 
rates the effect of confinement on both open and closed 
channels. Even in the degenerate, strongly interacting 
quantum gas on resonance, we observe a long lifetime, 
making our system promising for observing two-orbital 
Fermi gases with exotic order parameters [7]. Further¬ 
more, the tunability of the |e)-jg) interaction strength 
suggests novel avenues for the experimental implementa¬ 
tion of two-orbital many-body lattice models [36]. 
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TWO-CHANNEL MODEL IN A TRAP 


To model two atoms of mass m on a single site, we first consider the two-body problem in a 3D harmonic 

trap, V{r) = In this case, the center-of-mass and relative coordinates can be decoupled, and in the relative 

basis, we have Hamiltonian % = Hq -\- V, with 


Ho = '^en |o,n) (o,n| + ^ (e„ -b A^iB) \c,n) (c,n| (S.l) 

n n 

and interaction part 

V ^'^(Pn{0)(pn>{0) {u~^ \+,n) {-\-,n'\ -\-U~ \-,n) (S.2) 

n,n' 


Here, n labels the relative harmonic oscillator states with angular momentum I = 0 (these are the only ones that are 
affected by the short-range interactions), (pn(0) is the real-space harmonic oscillator wavefunction at r = 0, and the 
non-interacting energy e„ = 2nhujr- The triplet and singlet configurations are defined as |±) = ^ (|o) ± |c)), with 
corresponding interaction strengths U^, which do not depend directly on the nuclear spin. For simplicity, we focus 
on a zero-range interaction here, but finite-range corrections are easily incorporated into the formalism provided the 
range of the interactions is much smaller than the harmonic oscillator length = ^yh/mcor- 

To determine the two-body eigenstates, we consider the general wavefunction in the singlet-triplet basis 


1^) = 1 + ’ 1“’ ^)) 

n 


(S.3) 


and then insert this into the Schrodinger equation as follows: 


l(l l\(e^-E 0 \ 

2\1 -IJ \ 0 en + ApLB-Ej 




(S.4) 


By summing over n and replacing the bare interactions with the scattering lengths a^, we arrive at the matrix 
equation for the regular part of the two-channel wave function 


where 


and the quantity 
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contains the pair propagator in a harmonic potential [SI] 


U(E) = r{-E/2hwr) 

y/mlr T{—E/2hoJr — 1/2) 


(S.8) 


The interaction energies are thus obtained from the condition det (tq ^ + n(_E)) = 0. 

For finite-range corrections, we need to replace Tq ^ by an energy-dependent matrix that depends on the scattering 
phase shifts, as in Eq. (1) of the main text: 


r {E^)-^E^^ ^ cot6-{E^) 


(S.9) 


where the collision energy is Ec = E — We are interested in the leading order terms of the low-energy 

expansion: 


^/mEc 

n 


cot 6fgiE^) 



1 mEc 


(S.IO) 


where denotes the effective range. Formally, such an effective range may be included by using a two-channel model 
for each of the singlet and triplet interactions. For the trapped two-body problem, the low-energy expansion of r is 
thus equivalent to making the replacement l/afg >->■ l/afg — (^E — in Eq. (S.6). 

The interaction energies of the two-body problem in a 3D harmonic trap are dependent on the trap energy spectrum 
for the single-channel problem [Slj. In order to take into account the anharmonic character of the lattice potential, 
we expand this potential to the 6th order and correct the harmonic energy spectrum using second-order perturbation 
theory. This result matches very well with the diagonalization of the full Hamiltonian performed by Deuretzbacher 
et. al. [S2]. 

Using the determined scattering lengths and the corresponding effective ranges of the interactions, one can calculate 
the scattering amplitude resonance position in the bulk Bq = jebl/A/r = 42g^ G and the zero crossing Ba = 319g2 G. 
These values correspond to zero energy scattering and the error bars were obtained allowing 10% of uncertainty in 
the determination of the scattering lengths. 


EXPERIMENTAL METHODS 

In the cross-dimensional thermalization measurement, the atomic cloud is transferred into a crossed dipole trap 
operating at the magic wavelength with trap frequencies = Stt x (20,120,160) Hz. Simultaneously, the magnetic 
field is changed in a two-step ramp: First it reaches a value 100 G above the value of interest in 100 ms and then 
is finally ramped with a fixed speed of 2G/ms. Directional heating is produced in a controlled way by flashing the 
lattice potential along the vertical axis. Twenty pulses are applied over a time of three trap oscillation periods, in 
order to ensure that the cloud is thermalized along the excitation direction. After a variable hold time, the trap is 
switched off and the atoms in the ground state are imaged after 12 ms of free flight. 

In order to calibrate the magnetic field in our experiments, a measurement of the Zeeman shift of the 
I^SojTOi;’ = ±5/2) —)• pPi,m_F = ±7/2) transitions between stretched states is performed. A subsequent measure¬ 
ment of the Zeeman shift of the | ^Sq, tuf = ±5/2) —)■ pPo, rnp = ±5/2) transitions showed very good agreement with 
previously published values [S3, S4]. 


INTERORBITAL CROSS THERMALIZATION 

As described in the main text, a cross-dimensional thermalization experiment is performed to determine the 
magnetic-field dependence of the elastic scattering cross-section of a |g t) |e i) mixture. The aspect ratio of the 
cloud is measured by performing a Gaussian fit to the density distribution obtained by absorption imaging after 12 ms 
of time of flight. As an example, the time evolution of the cloud sizes for two magnetic field values is plotted in 
Fig. SI. The magnetic-field independent thermalization of the \g f) \g /) mixture is shown for comparison in Sl(a). 

For the cross-thermalization rate measurement, two different values of atomic density are used. For magnetic fields 
below 200 G, close to the Feshbach resonance, the density is lower in order to compensate for the large scattering cross- 
section. At higher magnetic fields, the higher value of atomic density is used. The two data sets are merged together 
using as reference the corresponding magnetic-field independent thermalization rates of the \g f) \g /) mixture. 
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FIG. SI. Time evolution of the cloud widths for two different magnetic fields for (a) a mixture of 1^ t) \g i) atoms and (b) a 
mixture of \g t) |e4,) atoms. Inset of (a): Trap geometry of thermalization experiment, with dipole trap beams in light red, 
heating lattice beams in dark red and the imaging beam in blue. The images are oriented in the 1 / 2 -plane. 


[51] T. Busch, B. Englert, K. Rziizewski, and M. Wilkens, Found. Phys. 28, 549 (1998). 

[52] F. Deuretzbacher, K. Plassmeier, D. Pfannkuche, F. Werner, C. Ospelkaus, S. Ospelkaus, K. Sengstock, and K. Bongs, 
Phys. Rev. A 77, 032726 (2008). 

[53] N. Poll, Z. W. Barber, N. D. Lemke, C. W. Oates, L. S. Ma, J. E. Stalnaker, T. M. Fortier, S. A. Diddams, L. Hollberg, 
J. C. Bergquist, A. Brusch, S. Jefferts, T. Heavner, and T. Parker, Phys. Rev. A 77, 050501 (2008). 

[54] F. Scazza, C. Hofrichter, M. Hofer, P. C. de Groot, I. Bloch, and S. Foiling, Nat. Phys. 10, 779 (2014). 









